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Abstract 

Quantum simulations of complex fermionic systems suffer from a variety of challenging problems. 
In an effort to circumvent these challenges, simpler "semi-classical" approaches have been used to 
mimic fermionic correlations through a fictitious "Pauli potential" . In this contribution we examine 
two issues. First, we address some of the inherent difficulties in a widely used version of the Pauli 
potential. Second, we refine such a potential in a manner consistent with the most basic properties 
of a cold Fermi gas, such as its momentum distribution and its two-body correlation function. 
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I. INTRODUCTION 



Insights into the complex and fascinating dynamics of Coulomb frustrated systems across 
a variety of disciplines are just starting to emerge (see, for example, Ref. [l| and references 
therein). In the particular case of neutron stars, one is interested in describing the equation 
of state of neutron-rich matter across an enormous density range using a single underlying 
theoretical model. In recent simulations we have resorted to a classical model that, while 
exceedingly simple, captures the essential physics of Coulomb frustration and nuclear sat- 
uration jH, 0, |4|. The model includes competing interactions consisting of a short-range 
nuclear attraction (adjusted to reproduce nuclear saturation) plus a long-range Coulomb 
repulsion. The charge-neutral system consists of electrons, protons, and neutrons, with the 
electrons (which at these densities are no longer bound) modeled as a degenerate free Fermi 
gas. 

So far, the only quantum effect that has been incorporated into this "semi-classical" 
model is the use of an effective temperature to simulate quantum zero-point motion. The 
main justification behind the classical character of the simulations is the heavy nature of the 
nuclear clusters. Indeed, at the low densities of the neutron-star crust, the de Broglie wave- 
length of the heavy clusters is significantly smaller than their average separation. However, 
this behavior ceases to be true in the transition region from the inner crust to the outer core. 
At the higher densities of the outer core, the heavy clusters are expected to "melt" into a 
collection of isolated nucleons with a de Broglie wavelength that becomes comparable to 
their average separation. Thus, fermionic correlations are expected to become important in 
the crust-to-core transition region. Unfortunately, in contrast to classical simulations that 
routinely include thousands — and even millions — of particles, full quantum-mechanical 
simulations of many-fermion systems suffer from innumerable challenges (see Ref. [g] and 
references therein). In an effort to "circumvent" - - although not solve — some of these 
formidable challenges, classical simulations of heavy-ion collisions and of the neutron-star 
crust have resorted to a fictitious "Pauli potential" . Within the realm of nuclear collisions, 
the first such simulations were those of Wilets and collaborators p]. Other simulations with 
a more refined Pauli potential have followed @, S @, ID, 11], but the spirit has remain 



the same: introduce a momentum dependent, two-body Pauli potential that penalizes the 
system whenever two identical nucleons get too close to each other in phase space. 

A goal of the present contribution is to show that the demands imposed by such a Pauli 
potential are too weak to reproduce some of the most basic properties of a zero-temperature 
(or cold) Fermi gas. Thus, we aim at refining such a potential in a manner that three 
fundamental properties of a cold Fermi gas be reproduced. These are: (i) the kinetic energy 
(as others have done before us), (ii) the momentum distribution, and (Hi) the two-nucleon 
correlation function. 

The manuscript has been organized as follows. In Sec. [TTJ some of the basic properties of a 
free Fermi gas are discussed. As our classical simulations must by necessity be carried out at 
finite temperature, a Sommerfeld expansion is used to compute low-temperature corrections 
to these observables. The section concludes with a review of the "standard" form of the 
Pauli potential, its flaws, and the measures that we take to overcome these flaws. In Sec. Illfl 
results for the kinetic energy, momentum distribution, and two-body correlation function 
are presented and contrasted against exact analytic results. Conclusions and possible future 
directions are presented in Sec. [TV] 
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II. FORMALISM 



The present section starts with a brief overview of some fundamental properties of a low- 
temperature Fermi gas [12|. Next, a fictitious Pauli potential is introduced and constrained 
to reproduce these properties via a purely classical simulation. In particular, special em- 
phasis is placed on the necessary modifications to the "standard" Pauli potential that are 
required to reproduce such fundamental properties. 



A. Free Fermi Gas 

The zero temperature Fermi gas is the simplest many-fermion system. Such a system dis- 
plays no correlations beyond those imposed by the Pauli exclusion principle and is described 
by the following free Hamiltonian: 



N 



H = y^. (i) 



2m 
1=1 

Here m is the mass of the fermion and N denotes the (large) number of particles in the 
system. As no interaction of any sort exists among the particles, the eigenstates of the system 
are given by a product of (single-particle) momentum eigenstates, suitably antisymmetrized 
to fulfill the constraints imposed by the Pauli principle. For simplicity, we assume that the 
fermions reside in a very large box of volume V = L 3 and that the momentum eigenstates 
satisfy periodic boundary conditions. We will be interested in the thermodynamic limit of 
iV —>■ oo and V —>■ oo, but with their ratio fixed at a specific value of the number density 
p = N/V. 

The ("box") normalized momentum eigenstates are simple plane waves. That is, 

V?p(r) = -^- r . (2) 



Given that the eigenvalue problem is solved in a finite box using periodic boundary condi- 
tions, the resulting single-particle momenta are quantized as follows: 

p(n) = — n = — {n x ,n y ,n z ) , with n { = 0, ±1, ±2, . . . (3) 

with the corresponding single-particle energies given by e(p) = p 2 /2m. 

Up to this point the spin/statistics of the particles has not come into play. We are 
now interested in describing the ground state of a system of N non-interacting, identical 
fermions and the resulting many-body correlations. Such a zero-temperature state is ob- 
tained by placing all particles in the lowest available momentum state, consistent with the 
Pauli exclusion principle. Using fermionic creation and annihilation operators satisfying the 



following anti-commutation relations 13 



[A p ,Al] = <W and {A p , A p ,} = {4, = , (4) 
the ground state of the system may be written as follows: 



Pf 



|$FG> = n A !>l $ vac>, (5) 
p=0 
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where |$ V ac) represents the (non- interacting) vacuum state and the Fermi momentum pp 
denotes the momentum of the last occupied single-particle state. Note that henceforth, no 
intrinsic quantum number (such as spin and/or isospin) will be considered. In essence, one 
assumes that all intrinsic degrees of freedom have been "frozen" , thereby concentrating on 
a single fermionic species (such as neutrons with spin up). In what follows, we compute 
expectation values of various quantities in the Fermi gas ground state (|$ vac )). 

We start by computing the Fermi momentum p F in terms of the number density of the 
system p = N/V. That is, 



N 



n FD (nj 



V 



d 3 p 



r^FDlPJ 



= V A 

=0 67T 2 



(6) 



or equivalently, 



Pf 



(6tt 2 p) 



1/3 



(7) 



Note that in Eq. ([6]) npo(p) denotes the Fermi-Dirac occupancy of the single-particle state 
denoted by p (or n) and the thermodynamic limit has been assumed. As all ground-state 
observables will be computed over a spherically-symmetric Fermi sphere, we define the Fermi- 
Dirac momentum distribution f(q) as follows: 



/(<?) = 3g 2 n FD (<?) , with / f(q)dq = 1 



(8) 



where the dimensionless quantity q=p/p F is the momentum of the particle in units of the 
Fermi momentum. 

All classical simulations performed and reported in the next sections must be carried out 
by necessity at finite temperature. Thus, we now incorporate finite temperature corrections 
to the various observables of interest. For temperatures T that are small relative to the 
Fermi temperature T F (with T F = ep), finite-temperature corrections may be implemented 
by means of a Sommerfeld expansion [l2| . For example, to lowest order in r = T/T F the 
momentum distribution becomes 



3q 2 



exp 



7T 



1 + — r 
12 



+ 1 



3g 2 6(l-g). 



(9) 



Here Q(x) is the "Heaviside step function" appropriate for a zero-temperature Fermi gas. 
Similarly, the energy-per-particle of a "cold" Fermi gas may be readily computed. One 
obtains [l2( 



E/N = e F J q 2 f(q)dq = -e F 



1 + 



5tt^ 
12" 



+ 0(r 4 



(10) 



where the Fermi energy is defined by €F = p 2 ? /2m. 

The last Fermi-gas observable that we focus on is the "two-body correlation function" 
g(r). This observable measures the probability of finding two particles at a fixed distance 
r from each other. Moreover, the two-body correlation function is a fundamental quantity 
whose Fourier transform yields the static structure factor, an observable that may be directly 
extracted from experiment. As such, the two-body correlation function is the natural meet- 
ing place of theory, experiment, and computer simulations 14|. The two-body correlation 
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function may be derived from the density-density correlation function 



131 



That is, 



9 r = — P2 x,y 
P 2 



-(^( x )^t (y) ^ (y) ^ (x) 



(11) 



where ip{x) is a fermionic field operator and the Dirac brackets denote a thermal expectation 
value. As the two-body correlation function for a non-interacting Fermi gas may be readily 
evaluated at zero temperature [13] , we only provide its extension to finite temperatures. To 
lowest order in r = T/Tp (and for a single fermionic species) one obtains 



g{r) = 1 



71 



2^2 



1 Z Z T 

12 



(12) 



where z=p^r and ji(z) is the spherical Bessel function of order n—1, namely, 



h{z) 



sin \z) 



cos z 



(13) 



Equations (jSJ), (|10|l . and (TP2]) display the three fundamental observables of a cold Fermi 
gas that we aim to reproduce in this work via a momentum-dependent, two-body Pauli 
potential. Note that in most (if not all) earlier studies of this kind, only the kinetic energy 
of the Fermi gas [Eq. ({TO]) ] was used to constrain the parameters of the Pauli potential p, 



a a a llOi , 111] . We are unaware of any earlier effort at including more sensitive Fermi- 
gas observables to constrain the parameters of the model. Clearly, it should be possible 
to reproduce the second moment of the distribution {i.e., the kinetic energy) even with 
an incorrect momentum distribution. Thus, while we build on earlier approaches, we also 
highlight some of their shortcomings. 



B. Pauli Potential: A New Functional Form 



In the previous section the wave function of a zero-temperature Fermi gas was introduced 
as follows: 



Pf 



(14) 



l*FG) = n A pl $ vac). 
p=0 

Essential to the dynamical behavior of the system are the anti-commutation relations 
[Eq. (j3j)] that enforce the Pauli exclusion principle {i.e., (AjJ 2 = 0). As it is often done, 
one may project the above "second-quantized" form of the many-fermion wave-function into 
configuration space to obtain the well-known Slater determinant. That is, 



$ FG (p 1 ,...,p N ;r 1 ,...,r N ) 



^Pll r l, 



Vpiir/vJ 



(15) 



where the single-particle wave- functions f p {r) are the ("box") normalized plane waves de- 
fined in Eq. (T5]). The Slater determinant embodies important correlations that were discussed 
in the previous section and that we aim to incorporate into our classical simulations. These 
are: 
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(a) At zero temperature, only momentum states having a magnitude |p| less than the Fermi 

momentum pp are occupied; the rest are empty. 

(b) As a consequence of the Pauli exclusion principle [(A^) 2 = 0], the probability that two 

fermions share the same identical momentum is equal to zero. Mathematically, this 
result follows from the fact that the wave-function vanishes whenever two rows of the 
Slater determinant are equal to each other. 

(c) Similarly, the wave-function also vanishes whenever two columns of the Slater determi- 

nant are identical. This fact precludes two fermions from occupying the same exact 
location in space. 

The first two properties are embedded in the momentum distribution of Eq. namely, a 
quadratic momentum distribution sharply peaked at the Fermi momentum p-p (recall that 
such a momentum distribution emerges after folding the Heaviside step function with the 
phase space factor). The third property induces spatial correlations that are captured by the 
two-body correlation function g(r) of Eq. ffl2l) . Indeed, the two-body correlation function 
is related to the integral of the Slater determinant over all but two of the coordinates of 
the particles (e.g., ri and r 2 ). Clearly, the Slater determinant vanishes whenever r x = r 2 , 
and so does the two-body correlation function at r = \ri — r 2 | = 0. It is the aim of this 
contribution to build a Pauli potential that incorporates these three fundamental properties 
of a free Fermi gas. 

However, before doing so, we briefly review the Pauli potential introduced by Wilets and 
collaborators — and used by others with minor modifications — to simulate the collisions 
of heavy ions and the properties of neutron rich matter at sub-saturation densities. Such a 
Pauli potential is given by a sum of momentum-dependent, two-body terms of the following 
form: 

JV 

Vpauii(Pi,---,Pjv;ri,...,r JV ) = V ° ex P(~4/ 2 ) ' ( 16 ) 

i<j=i 

where Vo > and the dimensionless phase-space "distance" between points (pi,i"j) and 
(Pj, r j) is g iven by 

Po 'o 

Here p$ and r are momentum and length scales related to the excluded phase-space volume 
that is used to mimic fermionic correlations. That is, whenever the phase-space distance 
between two particles is such that s?- < 1, then a penalty is levied on the system in an effort 
to mimic the Pauli exclusion principle. Although the parameters of this Pauli potential 
(Vo, po, and r ) can — and have — been adjusted to reproduce the kinetic energy of a free 
Fermi gas, it fails (as we show later) in reproducing more sensitive Fermi-gas observables, 
particularly, the momentum distribution f(p) and the two-body correlation function g(r). 

Upon closer examination, the above flaws should not come as a surprise. In our previous 
discussion of the Slater determinant it has been established that the probability of finding 
two identical fermions in the same location in space "or" with the same momenta must 
be identically equal to zero. Yet the Pauli potential of Eq. (1161) fails to incorporate this 
important dynamical behavior. Indeed, the above Pauli potential imposes a penalty on the 
system only when both the location "and" momenta of the two particles are close to each 
other (i.e., s^ <l). In particular, no penalty is imposed whenever two fermions occupy the 
same location in space, provided that their momenta are significantly different from each 



6 



other, i.e., |pj — Pj| 2 3>Po- Thus, the Pauli potential of Eq. ([161) will generate an incorrect 
two-body correlation function, namely, one with g(r) ^0 as r tends to zero. By the same 
token, an incorrect momentum distribution will be generated, although not necessarily its 
second moment. 

To remedy these deficiencies, a Pauli potential is now constructed so that the three 
properties [(a), (b), and (c)] defined above are explicitly satisfied. To this end, we introduce 
the following form for the Pauli potential: 

JV 



Vpauii(Pi,---,Piv;ri,...,r i v) = ^2 [v r A exp(-r ii /r ) + V B exp(-p i:j /p ) 

i<j=i 

N 

+ J2 V c®v(<li) > (18) 

2=1 

where — |r» — pij = \pi — pj\, Qi = \Pi\/pF, and Q v is a suitably smeared Heavisi de-step 
function of the following form: 

= TT r 1 ( 2 Tvi — (?) • ( 19 ) 

1 + exp| — n(q 2 — 1)\ v^co 

The parameters of the model Va,Vb,Vc and r ,p ,r] will be adjusted to reproduce both 
the momentum distribution and two-body correlation function of a low-temperature Fermi 
gas. Note that the phase-space dependence of the Pauli potential has been separated into 
a "sum" of two-body pieces, with the first acting exclusively in configuration space and the 
second one only in momentum space. The first term in the potential imposes a penalty as 
the particles get too close (of the order of r ) to each other. Similarly, the second term in 
the potential penalizes particles whenever their relative momenta becomes of the order of 
Pq. Finally, the third "one-body" term enforces the low-temperature behavior of the Fermi 
gas, namely, that the probability of finding any particle with a momentum significantly 
larger than the Fermi momentum is vanishingly small. Most of the parameters will depend 
explicitly on the density of the system (see Sec. IIIip . This reflects the complex many-body 
nature of the Pauli correlations and our inability to simulate them by means of a "simple" 
(albeit momentum dependent) two-body potential. 



III. RESULTS 



We start this section by listing in Tabled the parameters of the Pauli potential. For the 
strength parameters, the following simple density dependence is assumed: 

Vi{ P ) = V?(p/ PQ ) a \ z = {A,B,C} } (20) 

where po = psat/^ = 0.037 fm -3 is the density of a single fermionic species (for example, 
neutrons with spin up) at nuclear-matter saturation density (p sa t = 0.148 fm -3 ). For the 
range parameters (r and po) the following scaling relation is adopted: 

ro = Pa/pf , (21a) 
Po = PbPf ■ (21b) 
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V° 

V B 


V° 
V C 


OLA 


OLB 


OLC 


/3a 


Pb 


7/ 


13.517 


1.260 


3.560 


0.629 


0.665 


0.831 


0.845 


0.193 


30 



TABLE I: Strength (in MeV) and range parameters (dimensionless) for the various components 
of the Pauli potential. See Eqs. (fl~8|) . ([20]) . and (f2Tj) . These values have been used to simulate a 
system of A = 1000 identical fermions at a temperature of r = T/Tp = 0.05. 



Monte-Carlo simulations for a system of N = 1000 identical fermions at the finite (but 
small) temperature of r = T/Tp = 0.05 have been performed. Initially, the particles are 
distributed randomly throughout the box with momenta that are uniformly distributed up to 
a maximum momentum of the order of the Fermi momentum. After an initial thermalization 
phase of typically 2000 sweeps (or 2 million Monte Carlo moves for both coordinates and 
momenta) data is accumulated for an additional 2000 sweeps, with the data divided into 10 
groups to avoid correlations among the data. It is from these 10 groups that averages and 
errors are generated. 

A. Kinetic Energy 

The kinetic energy of a finite-temperature Fermi gas as a function of density is displayed 
in FigJTJ The (black) solid line represents the analytic behavior of a free Fermi gas correct 
to second order in r. The (red) line with circles is the result of the Monte Carlo simulations 
with the Pauli potential defined in Eq. ffT8l) . The agreement (to better than 5%) is as good 
as the one obtained with earlier parametrizations of the Pauli potential. 

To our knowledge, reproducing the kinetic energy of a free Fermi gas is the sole constraint 
that has been imposed on most of the Pauli potentials available to date. However, while a 
host of Pauli potentials can reproduce such a behavior, it is unclear if these potentials can 
also reproduce the full momentum distribution. Thus, we now show how the Pauli potential 
defined in Eq. (1181) is successful at reproducing two highly sensitive Fermi-gas observables, 
namely, the momentum distribution and the two-body correlation function. 

B. Momentum Distribution 

The momentum distribution obtained from the Monte-Carlo simulations is displayed in 
Fig. [2] for a variety of densities. Note that the momentum distribution has been normalized 
to one and that the densities have been expressed in units of po = 0.037 fm -3 . As indicated 
in Eq. (Q, the momentum distribution of a cold Fermi gas depends solely on the two di- 
mensionless ratios q = p/pF and r = T/Tp. Thus, all curves must collapse into the exact one 
- displayed by the (black) solid line — independent of density. It is gratifying to see that 
this is indeed the case. In contrast, we show in the next section that the standard Pauli 
potential of Eq. (TTBT) fails to reproduce this behavior. 
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FIG. 1: (color online) Average kinetic energy of a system of iV = 1000 identical fermions at a 
temperature of r = T/Tp = 0.05. The line with circles is the result of the Monte Carlo simulations 
with the Pauli potential of Eq. (|18p . The solid line is the exact behavior of a non-relativistic Fermi 
gas, as given by Eq. (fTUj) , 



C. Two-Body Correlation Function 

The two-body correlation function of a free Fermi gas is displayed in Fig. [3) The two-body 
correlation function g(r) is related to the probability of finding two particles separated by 
a distance r. For a free Fermi gas, the probability of finding two identical fermions at zero 
separation is identically equal to zero [see Eq. ( |T2l) ]. As this short-range (anti-) correlation 
is the sole consequence of the Pauli exclusion principle, the "Pauli hole" disappears for 
distances of the order of the interparticle separation (pp 1 ). As in the case of the momentum 
distribution, the correlation function depends only on two dimensionless variables (z=ppr 
and r). It is again gratifying that the simulation curves scale to the exact correlation 
function, depicted here with a (black) solid line. 



D. Comparison to other approaches 

In this section we compare the Pauli potential introduced in Eq. (fT8|) to earlier approaches 
that are based on Eq. (jl~6l) . Such approaches have been very successful in reproducing the 
kinetic energy of a free Fermi gas for a wide range of densities. Indeed, the kinetic energy 
displayed in Fig. [2] of Ref. [l(| is as good — if not better — than the one obtained here. 
However, a faithful reproduction of the kinetic energy does not guarantee that the system 
displays the same phase-space correlations as that of a free Fermi gas. To illustrate this point 
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FIG. 2: (color online) Momentum distribution of a system oi N = 1000 identical fermions at a 
temperature of r = T/Tp = 0.05 for a variety of densities (expressed in units of po = 0.037 fm~ 3 ). 
The momentum distribution has been normalized to one [see Eq. (jHJ)]. The (black) solid line with 
no symbols gives the exact behavior of a non-relativistic Fermi gas. 



we compare in Fig. H] the results from the two approaches for the momentum distribution 
and two-body correlation function at a fixed density of po = 0.037 fm~ 3 . The left-hand panel 
in the figure displays the momentum distribution and indicates that while earlier approaches 
(depicted with a blue line with circles) are accurate at reproducing the second moment of the 
distribution (i.e., the kinetic energy) the momentum distribution itself shows a behavior that 
differs significantly from that of a cold Fermi gas. The right-hand panel shows deficiencies 
that are as — or even more — severe. Two points are worth highlighting. First, the two-body 
correlation function g(r) generated with the standard form of the Pauli potential does not 
vanish at r = 0. Second, for distances of the order of the inter-particle separation and beyond, 
the two-body correlation function develops artificial oscillations. Failure in reproducing the 
correct behavior of g(r) at r = is relatively simple to understand. Potentials based on 
Eq. ffTBT) impose a penalty on the system only if both the positions and momenta of the 
two fermions are close to each other. Yet the correlations embodied in a free Fermi gas 
are significantly more stringent than that. Indeed, a Slater determinant vanishes if either 
the positions or the momenta of the two fermions are equal to each other. In contrast, the 
development of artificial structure in g(r) is a more subtle effect that is intimately related 
to the momentum dependence of the potential and that we now address. 

The artificial oscillations present in g(r) are reminiscent of the structure of liquids and/or 
crystals. The emergence of crystalline structure in a low-temperature/low-density system 
would be expected if the energy of localization becomes small relative to the mutual repulsion 
between the particles. Such, however, is not the behavior of a cold Fermi gas. While spatial 
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z=p F r 



FIG. 3: (color online) Two-body correlation function for a system of iV = i000 identical fermions at 
a temperature of r = T/Tp = 0.05 for a variety of densities (expressed in units of po = 0.037 fm~ 3 ). 
The (black) solid line with no symbols gives the exact behavior of a non-relativistic Fermi gas. 



anti-correlations would favor the formation of a periodic structure, the large velocities of 
the particles — resulting from the Pauli exclusion principle — would make their localization 
extremely costly. What is not evident is the reason for the momentum distribution generated 
with a standard Pauli potential to lead to crystallization, whereas not that of a real Fermi 
gas (see left-hand panel in Fig. H]). 



To our knowledge, the answer to this question was first provided by Neumann and Fai [15 



One must realize that any Hamiltonian that contains momentum-dependent interactions - 
such as most (if not all) Pauli potentials — generates a ''canonical" momentum distribution 
that may (and does!) differ significantly from the corresponding "kinematical" momentum 
distribution (tt) = (mr). Indeed, in a Hamiltonian formalism where the Hamiltonian de- 
pends on the positions and canonical (not kinematical) momenta of all the particles, namely, 
H = H(pi, . . . , pjv; ri, . . . , rjy) the kinematical velocities must be obtained from Hamilton's 
equations of motion, i.e., = dH/dpi. For a Hamiltonian that contains momentum- 
dependent interactions — as in the case of the Pauli potential — then the kinematical 
momentum r K i = mi i differs from the corresponding canonical momentum p^. This suggests 
that while a given choice of Pauli potential may produce the correct canonical momentum 
distribution, it may generate kinematical velocities that are too small to prevent crystal- 
lization. Such a distinction between canonical and kinematical momenta is an unwelcome, 
yet unavoidable, consequence of the approach. A cold Fermi gas — a system subjected to 
no interactions — is the quintessential quantum system were such an artificial distinction is 
not required. 

It is worth noting that although some earlier choices yield an extremely soft kinematical 
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FIG. 4: (color online) Comparison between the Pauli potential introduced in this work [Eq. (|18p ] 
and earlier approaches based on Eq. (|16p . The left-hand panel shows the momentum distribution 
while the right-hand panel the two-body correlation function. The simulations have been carried 
out for a system of N = 1000 identical fermions at a density of po = 0.037 fm~ 3 . The (blue) line with 
circles is obtained from a Monte Carlo simulation using Eq. (|16p with the parameters of Ref. [10]. 
The (red) line with squares displays the results using the Pauli potential introduced in this work; 
the black line gives the exact behavior of a non-relativistic Fermi gas. 



momentum distribution, the Pauli potential introduced in this work is not immune to such a 
disease. Indeed, we believe that an unrealistically soft kinematical momentum distribution 
is likely to be a general result of the Hamiltonian approach. Thus, we close this section by 
displaying in Fig. the kinematical velocity distribution obtained with the new choice of 
Pauli potential introduced in Eq. (fTSlb For comparison, the exact distribution of a cold Fermi 
gas (solid black line) is also included. The dependence of the Pauli potential on the canonical 
momentum distribution is responsible for generating such a soft velocity distribution, with 
its peak around 1/10 of the Fermi velocity. And while we were able to avoid crystallization 
with the present set of parameters (see Figs. [3] and H]), the risk of crystallization looms large 
(see next section). 



E. Finite-Size Effects 

We have observed that the Pauli potential introduced in Eq. (fl~8{) . with its parameters 
suitably adjusted, has been successful in reproducing a variety of Fermi-gas observables, such 
as its kinetic energy, its (canonical) momentum distribution, and its two-body correlation 
function. However, the first indication of a potential problem — and one that may be 
generic to all approaches employing momentum- dependent interactions — is the emergence 
of an unrealistically soft velocity distribution and with it, the possibility of artificial spatial 



12 



9.0 



-• p=0.2 
-■ p=0.4 
p=0.6 
^p=0.8 
-rp=1.0 



N=1000 
T =T/T =0.05 

F 




u=v/v 



FIG. 5: (color online) "Kinematical" velocity distribution of a system of N = 1000 identical fermions 
at a temperature of r = T/Tp = 0.05 for a variety of densities (expressed in units of po = 0.037 fm~ 3 ). 
The velocity distribution has been normalized to one. The (black) solid line with no symbols gives 
the exact behavior of a non-relativistic Fermi gas. 



correlations (i.e., crystallization). Fortunately, with the choice of parameters adopted in 
this work, the problem of crystallization was avoided (see Fig. [3]). Yet, there is no guarantee 
that crystallization will not become a problem as one examines the sensitivity of our results 
to Unite-size effects. 

In order to estimate the sensitivity of our results to finite-size effects, Monte Carlo simu- 
lations were performed for a system containing iV = 250, iV = 500, iV=1000, and iV = 2000, 
identical fermions (note that the results reported so far have been limited to 1000 parti- 
cles). The conclusions from this study are mixed. First (and fortunately) no evidence of 
crystallization or of significant finite-size effects were found. These findings are displayed in 
Fig. [H] for both the canonical momentum distribution (left-hand panel) and the two-body 
correlation function (right-hand) panel. Unfortunately, however, in order to preserve the 
high quality of the results previously obtained with 1000 particles, a parameter of the Pauli 
potential [Vg in Eq. (|T8|) ] had to be fine tuned. Specifically, the following scaling with 
particle number was used: 

Vl(N) = V° B (N = im (^j , (22) 

with Vg(N = 1000) = 1.26 MeV being the value listed in Table [H This unpalatable fact may 
be a reflection of the highly challenging task at hand: how to simulate fermionic many-body 
correlations by means of a "simple" two-body Pauli potential. 
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FIG. 6: (color online) Finite-size effects on the canonical momentum distribution (left-hand panel) 
and the two-body correlation function (right-hand panel) for a system of identical fermions at a 
temperature of r = T/Tp = 0.05 and a density of p = po = 0.037 fm -3 . Simulations were carried 
out for systems containing N = 250, 500, 1000 and 2000 particles. The (black) solid line with no 
symbols gives the exact behavior of a non-relativistic Fermi gas. 



IV. CONCLUSIONS 



Semi-classical simulations of fermionic systems attempt to circumvent the many chal- 
lenges posed by bona-fide quantum simulations through the inclusion of a fictitious Pauli 
potential @, 0, I, I, M, ■ In the present contribution we conducted a critical study 
of the features and predictions of a widely used version of the Pauli potential. We con- 
cluded that the constraints imposed by such a Pauli potential, namely, the suppression of 
phase-space configurations having two fermions with both positions and momenta similar to 
each other, are too weak to faithfully reproduce some basic properties of a free Fermi gas. 
Instead, by examining the well-known behavior of the Slater determinant we suggest that 
phase-space configurations should be suppressed when either the positions or the momenta 
of the fermions are close to each other. By incorporating these features into a new form of 
the Pauli potential — and by carefully tuning the parameters of the model — the momen- 
tum distribution and the two-body correlation function of a free Fermi gas were accurately 
reproduced. 

In the course of this study a pathology that is generic to all Pauli potential was un- 
covered. As pointed out by Neumann and Fai 15| (to our knowledge for the first time), 
the "canonical" momentum distribution generated via Monte Carlo (or other) methods may 
differ significantly from the resulting "kinematical" momentum distribution. This suggests 
that while the kinetic energy of the free Fermi gas (computed from the canonical momenta) 
may be accurately reproduced, the distribution of velocities may be grossly distorted. In- 
deed, we found a distribution of velocities that significantly under-estimates — by a factor 
of 10 — that of a free Fermi gas. Such "sluggishness" among the particles could have dis- 
astrous consequences by inducing artificial ordering in the system (e.g., "crystallization"). 
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The possible appearance of artificial long-range order in the system must be examined on 
a case by case basis. For example, with the standard version of the Pauli potential [l(| the 
system displays an anomalous two-body correlation function suggestive of crystallization. 
On the other hand, the Pauli potential introduced in this work faithfully reproduces the 
two-body correlation function of a free Fermi gas. 

In conclusion a systematic study of the standard version of the Pauli potential has been 
conducted. While simple and widely used, such a version fails to reproduce some of the 
most basic properties of a free Fermi gas. Moreover, an analysis of the two-body correlation 
function generated by such a model reveals artificial long-range order. To correct these flaws 
a refined version of the Pauli potential was introduced. This new version — inspired by the 
properties of a Slater determinant — generated accurate (canonical) momentum distribution 
and two-body correlation functions, while avoiding crystallization. Yet the u kinematical" 
momentum distribution generated from such a Pauli potential was grossly underestimated. 
We believe this to be a generic feature of any momentum-dependent Pauli potential that is 
used in conjunction with a Hamiltonian approach Thus, the possibility for generating 
artificial correlations in the system {e.g., crystallization) looms large. 
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